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What is claimed is: 

1 . A method for determining whether any of a predetermined set of tones 
5 present in a plurality of successive frames of digital samples of a received signal 
falls within a predetermined frequency tolerance, comprising the steps of: 

obtaining discrete Fourier transform pairs of in-phase and quadrature dot 
products of said samples and integer multiples of a base frequency; 

determining the phase angle for the highest power ones of said products 
10 obtained on successive frames; and 

subtracting an offset phase from the phase angle of said highest power 
ones of said products to determine the deviation of an observed tone from said 
predetermined frequency tolerance. 

15 2. A method for determining said tones according to claim 1, wherein phase 
angles are computed by finding the arc tangent of said quadrature and in-phase 
dot products on successive ones of said frames. 

3. A method for determining said tones according to claim 2, wherein the 
20 phase angle corresponding to said arc tangent is approximated as the quotient of 

the quadrature and in-phase products for small absolute values of the quotient. 

4. A method for determining said tones according to claim 2, wherein said 
offset phase is proportional to the difference between an integer multiple of said 

25 base frequency and the closest one of said set of tones. 

5. A method for determining said tones according to claim 2, wherein a 
change in said phase angle between successive frames of said waveform is 
computed. 
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6. A method for determining said tones according to claim 5, wherein the 
amplitude of a detected tone is compensated by a multiplicative correction of 
x/Sin(x) based upon the frequency deviation (x) of the observed tone from a tone 
of said set. 
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7. A method for determining said tones according to claim 3, wherein said 
approximation of said arc tangent function Theta is performed as follows: 

a. for I >= 0, Q >= 0, Al > AQ, and 0 <= Theta <= 71/4, 

Theta = Thetal ; 

5 b. for I >= 0, Q >= 0, Al < AQ, and 71/4 <= Theta <= 7C/2, 

Theta = 71/2 -Thetal; 

c. for I < 0, Q >= 0, Al > AQ, and 3/471 <= Theta <= 71, 
Theta = 7C - Thetal ; 

d. for I < 0, Q >= 0, Al < AQ , and 7E/2 <= Theta <= 3/471, 
10 Theta = 71/2 + Thetal ; 

e. for I >= 0, Q < 0, Al > AQ, and -7T/4 <= Theta <= 0, 
Theta = -Thetal; 

f. for I >= 0, Q < = Al < AQ, and -7C/2 <= Theta <= -71/4, 

Theta = -7C/2 + Theta; 
15 g. for l< 0, Q < 0, Al > AQ, and -7T < Theta <= -3/471, 

Theta = -n + Thetal ; and 
h. for I < 0, Q < 0, Al < AQ, and -3/47E <= Theta <= -7C/2, 

Theta = -7C/2 - Thetal ; 

where I is the in-phase component; Q is the quadrature component; Al is 
20 the absolute value of I; AQ is the absolute value of Q; Theta 1 is the absolute 
value of AQ/AI. 



